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ABSTRACT 

Motivated by recent evidence that the epoch of reionization of hydrogen may 
have ended at a redshift as low as z ~ 6, we consider the detectability of the 
sources responsible for this reionization. The main idea is that reionization places 
limits on the mean surface brightness of the population of reionization sources. 
Reducing the number of model-dependent assumptions to a minimum, we dis- 
cuss the observability of these sources with existing and planned telescopes. We 
define a family of models characterized by two parameters: the Lyman contin- 
uum escape fraction from the sources, and the dumpiness parameter C of 
the intergalactic medium. The minimum surface brightness model corresponds 
to a value of unity for both parameters. We find that the detection of the non- 
ionizing UV continuum of the reionization sources will be difficult to accomplish 
before the launch of JWST if these sources have a mean surface brightness close 
to the minimum value. However, if the values of fc and C are more realistic, 
the reionization sources may well be detected by HST before JWST, perhaps by 
exploiting gravitational lensing amplification by foreground clusters of galaxies. 
Instead of a detection in the continuum, one may attempt to detect the Lya 
emission line by narrow-band imaging. Present Lya searches aX z > Q suggest 
that either the typical sources are fainter than dwarf galaxies, or that the escape 
fraction of ionizing photons is much higher than 50%, so that there are no bright 
compact HII regions formed around the ionizing sources. 

Subject headings: cosmology: reionization; high-z galaxies 



1. Introduction 

The recent discovery of quasars at 2; > 6 in the Sloan Digital Sky Survey (Fan et al. 
2001) has given a new impetus to the study of the epoch of reionization. In particular, SDSSp 
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J103027.10 +055255.0 ai z — 6.28 shows a significant continuum absorption shortwards of 
Lya, reminiscent of a Gunn-Peterson (1965) trough (Becker et al. 2001; Fan et al. 2002), 
and suggests that we may be close to identifying - or may have already identified - the end of 
the epoch of hydrogen reionization. A high priority now is to identify the sources responsible 
for producing the UV ionizing radiation (the reionization sources) and to determine their 
nature. There have been several attempts to predict the properties of these sources (Madau 
et al. 1999; Barkana and Loeb 2000; Miralda-Escude et al. 2000; Haiman et al. 2001; Oh 
et al. 2001; Bromm and Clarke 2002; Yan et al. 2002; Cen 2003; Pukugita and Kawasaki 
2003; Venkatesan et al. 2003; Wyithe and Loeb 2003a,b) and several reviews have been 
pubhshed recently (Loeb and Barkana 2002; Stiavelli 2002). 

Here, we introduce a new approach to explore the observable properties of the reion- 
ization sources, while minimizing the number of model-dependent assumptions. Our main 
interest is to determine the minimum cumulative surface brightness (in flux, maximum in 
magnitudes) of the population of reionization sources rather than to provide a complete 
evolutionary model of these sources. A strict lower limit to the surface brightness of the 
sources is obtained by requiring that they produce as many ionizing photons as there are 
hydrogen atoms in the universe. This limit can be improved upon by considering the main 
physical effects taking place during reionization and their relative importance. The required 
surface brightness is minimum when the escape of ionizing UV radiation is the highest, so 
that the source luminosity is lowest in the non-ionizing UV and in the visible. A maximum 
surface brightness is obtained by requiring that the reionization sources do not overproduce 
heavy elements. Our general approach is applicable to most types of reionization sources, 
but when specific numerical examples are given, we illustrate our procedure by focussing on 
population III stars, as these objects have very high effective temperatures and therefore are 
very effective producers of ionizing UV photons. Our mean surface brightness estimates are 
compared to the parameter space that can be probed by existing and future telescopes, in 
order to help planning the most effective surveys. 

In Section 2, we derive the minimum mean surface brightness of the reionization sources 
by maximizing the production and the release of Lyman continuum photons into the inter- 
galactic medium (IGM). We first recover the proportionality between surface brightness of 
the reionization sources and the comoving density of hydrogen (2.1). The proportionaUty 
constant depends on both the helium fraction and the spectral energy distribution (hereafter 
SED) of the reionization sources. We start out with the case of a pure-hydrogen IGM (2.2) 
and then consider the effects of primordial helium (2.3). The least model-dependent results 
in this paper are those for the homogeneous case (2.4), and the corresponding minimum 
metallicity enrichment of the universe (2.5). In Section 3, we consider more realistic models 
in order to provide a better comparison with observations. In particular, we discuss the 
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effects of a limited escape fraction of ionizing UV photons (3.1), increased recombinations 
induced by density inhomogeneities in the IGM (3.2), a luminosity function for the reioniza- 
tion sources (3.3), and the detectability of Lya emission from these sources (3.4). Finally, 
in Section 4, we discuss the observability of these sources both from space, either with HST 
or JWST, and from ground-based facilities. 



2. The minimum surface brightness for reionization 

2.1. The surface brightness of reionization sources 

There is a simple connection between the comoving density of hydrogen and the mean 
surface brightness Jy of the reionization sources in their non-ionizing continua. Neglecting 
the absorption of non-ionizing photons in the IGM, the surface brightness of a class of sources 
is simply related to their volume emissivity E{v), the power radiated per unit frequency per 
unit comoving volume, by 



J^^-^j^ E{[l + z{t)]u)dt, (1) 

where ti indicates the cosmic time when the sources turn on and t2 the time when reionization 
is complete (see, e.g., Peebles 1993). Here and elsewhere we refer to comoving volumes and 
densities unless otherwise specified. The production rate of ionizing photons per unit volume 
fic is related to the volume emissivity by 



where h is Planck's constant and vh = 2.467 x 10^^ Hz is the frequency of photons at the 
hydrogen ionization threshold. Here and in the following we denote with the subscript c 
quantities referring to Lyman continuum photons. 

For a given SED of reionization sources, there is a direct proportionality relation be- 
tween the production rate of ionizing photons ric and the emissivity E{i>q) at some reference 
frequency Ug. With this in mind, we define the proportionality coefficient 

AM^^. (3) 

Tic 

In the following, we choose uq to be a fixed frequency in the observer frame in the non-ionizing 
continuum of the sources. 
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In order to ionize hydrogen completely, the total number of ionizing photons per co- 
moving volume must exceed the comoving number density of hydrogen uh by some factor 
B that depends on the details of the recombination process and on the helium fraction 

1 

B = — / hcdt. (4) 
Combining equations (1), (3), and (4), we obtain 

cA{v)B 

Jv = — Uh, (5) 

where A{v) denotes the average of A{i'q) weighted by the production rate of ionizing photons, 
namely 



!l:A{ [l-rz{t)]v )hJt 
Kl ^cdt 



Equation (5) establishes a simple, direct relation, between the mean surface brightness of 
reionization sources and the comoving density of hydrogen atoms. In the following, for 
simplicity, we will neglect possible variations in fic- In this case, Alu) becomes a time 
average: 

A(u)^-^ rA{[l + z{t)]u)dt. (7) 

The volume emissivity E{i>) for a population of sources is given in terms of their lumi- 
nosity function per comoving volume, ^{L^), as 



POO 

E{v) = / ^L,)L,dL,. 
Jo 



Equations (1) and (8) establish the connection between the mean surface brightness of 
a population of sources and their luminosity function. In the case of identical sources the 
luminosity function is simply a delta-function. The following subsections will be devoted to 
estimating the values oi A{u) and B. 
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2.2. Reionization in a hydrogen-only IGM 

Due to recombinations, the minimum number of ionizing photons per comoving volume 
needed to reionize hydrogen exceeds the number of hydrogen atoms {B > 1). The equation 
for the ionized fraction x = Up/nn, including recombinations is 

^^i!^-CaBx{l + z{t)fn,, (9) 
at Uh 

where Ue is the comoving density of electrons, C =< rip > / < Up >^—< nl > / < Ue >^ 
is the IGM dumpiness factor, and as = 1.4 x 10~^^ cm^ s~^ is the recombination rate for 
an electron temperature oi T^, — 2 x 10*^ K (Pengelly 1964). This value of the electron 
temperature is appropriate for primordial HIT regions (Daub 1963; Davidson and Kinman 
1985; Tumlinson et al. 2001). The factor (1 + z)"^ is due to the fact that recombinations 
depend on the proper (not comoving) density of electrons. 

For the pure-hydrogen composition adopted in this section we have Ue = Up. Solving 
equation (9) and requiring x{t2) = 1, we can determine function of ti and ^2- 

The parameter B of equation (4) is the ratio of the actual number density of ionizing 
photons to the number density of hydrogen atoms. We define Bh as the value of B in the 
pure hydrogen case. In the absence of recombinations (i.e. in the limit of a very low density) 
and in the pure hydrogen case, the required comoving number density of ionizing photons, 
^c,noHe = ^c(^2 — ti) , equals the comoving number density of hydrogen, nc,noHe = and 
therefore Bh — 1- When recombinations are considered in the homogeneous case {C = 1), we 
find Bh < 1-03 or < 1.08 for Z2 — z{t2) ^ 10, and Az = z{ti) — z{t2) = 1 or 3, respectively. 
We also find that if reionization occurs at Z2 — 6, Bh < 1-3 for Az < 50. Thus, in a 
homogeneous (C = 1) pure-hydrogen IGM, recombinations are only a minor correction. 



2.3. Reionization in a hydrogen- helium IGM 

More reahstically, the chemical composition of the universe at reionization includes both 
hydrogen and hehum, with mass fractions X — 0.76 and Y — 0.24 (Pagel 2000). Thus, the 
mean molecular weight is given by {mHY / mHe-\- X)~^mH — l-^niH, where rriH and niHe are 
the masses of hydrogen and hehum atoms, respectively. Assuming = 0.268, Qb — 0.044, 
Qj^ = 0.732, and Hq = 71 km s"^ Mpc"^ (the best-fit WMAP parameters; Spergel et al. 
2003), we find that the number density of hydrogen atoms (including ions) per comoving 
volume, nn is 2.1 x 10~'^cm~^ or 6.2 x lO^^Mpc"^. 

Photons with energies above 24.4 eV can ionize He as well as H. The cross section for 
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He ionization at the threshold is higher than that of H by about a factor 7 but the number 
density of He is lower by a factor of about 12, so that photons at these energies may be 
absorbed by, and ionize, either atoms. In the following we will assume that He absorbs as 
many photons as possible in proportion to its abundance relative to H, i.e. about 8 % of 
the total budget of ionizing photons. Photons with energy above 54.4 eV can ionize helium 
twice. Whether helium is ionized or not and how many times, depends on the hardness of the 
UV spectrum. As an example, a blackbody at T = 10^ K emits 60% of its ionizing photons 
at energies greater than 24.4 eV. In this case, all hehum will be ionized, and about 10 % 
will be twice ionized. Therefore, we have Ue — l.OQup, and the required minimum number 
of ionizing photons will be higher by 9% relative to the pure H case. Note that if the SED 
were hard enough to double-ionize all helium, would increase to 1.16np. This is unlikely 
to happen for stellar sources, but it might for active galactic nuclei (AGNs). 

Helium ionization has an additional, indirect effect on the required number of ioniz- 
ing photons because the 9% increase in the electron number density makes recombinations 
slightly more effective. This can be estimated by integrating equation (9) for rig — 1.09np, 
from which we obtain 

B ~ 1.09[1 + 1.09{Bh - 1)]Bh. (10) 

In general, we find that primordial helium recombinations increase B by less than an 
additional 1% for Z2 < 10 and Az < 3 and by less than 2% for reionization at Z2 = 6 and 
Az < 50. Thus, recombinations remain a minor correction in the homogeneous (C = 1) 
case, even for a hydrogen- helium IGM. 

2.4. Results for a homogeneous IGM 

So far, the only dependence of our calculations on the type of ionizing sources has 
been in the assumption that all helium is ionized. Indeed, even if this assumption were not 
verified, our results would change only by 10 % or less. In order to derive the mean surface 
brightness of sources from their production of ionizing photons, we need to assume a specific 
SED and compute the quantity A{i/). For definiteness, we will now consider population III 
stars as the sources responsible for reionization of hydrogen. Pop III stars are characterized 
by virtually zero metallicities and consequently very high effective temperatures (Castellani 
et al. 1983; Tumhnson and Shull 2000; Marigo et al. 2001; Bromm et al. 2002). Here 
we will describe them as blackbodies with effective temperatures of 10^ K. Note that high 
temperatures around 10^ K are reached even if the mass function does not extend beyond 
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lOOM© (Tumlinson and ShuU 2000; Baraffe et al. 2001; Schaerer 2002). Pop III stars thus 
provide us with the minimum mean surface brightness needed for the case of reionization 
by stars. In fact, for a fixed flux of ionizing photons, cooler star arc brighter in the visible 
and in the non-ionizing UV continuum. One can define the efficiency of an ionization source 
as proportional to the ratio of the number of ionizing photons to the total energy output. 
Cooler stars emit a smaller fraction of their radiation in the form of ionizing photons and are 
therefore less efficient ionizers. As an example, for equal ionizing fluxes, a stellar population 
with = 5 X 10^ K is a factor ~4.5 brighter at A = 1400 A than a stellar population with 
= 10^ K. A similar comparison can be made for typical AGNs. We find that for typical 
power-law SEDs, AGNs are less efficient ionizers than Pop III stars. 

Once the SED is specified, we can determine the value of A from equation (3) and (7). 
In integrating equation (7), the emissivity can be assumed to be zero below the rest-frame 
wavelength of 1216 A since during the reionization era these photons will be absorbed by 
intervening neutral hydrogen. In magnitude units, the surface brightness per square arcmin 
is given by /iab — —2.5logJi, + 8.9, when Ji, is expressed in Jy arcmin"^. 

In Table 1, we give the minimum surface brightness for reionization sources for a variety 
of terminal redshifts of reionization Z2 and durations Az. Two values of surface brightness 
are given in the table for each combination of Z2 and Az, corresponding to the rest-frame 
1400 A wavelength observed at the terminal redshift of reionization Z2 or at the mean redshift 
z = {zi + Z2)/2. Thus, they arc computed at the observer wavelengths A = 1400(1 + Z2) A 
and A = 1400(1 + ^) A , respectively. In addition, we also give the number of photons per 
ionization B and the Thomson opacity tt, integrated to the beginning of the reionization 
process, and given by 



where ctt = 6.653 x 10 cm^ is the Thomson cross-section. 

All results in Table 1 include the effects of helium ionization (by population III stars) 
and recombinations. It is clear that in the homogeneous case (C = 1), the minimum sur- 
face brightness (in flux, maximum in magnitude) is essentially independent of the terminal 
redshift of reionization (within 0.1 magnitudes for Az = 1) but does depend on Az. 

The value of the Thompson opacity is a strong function of Az. In the case of extended 
reionization, Az = 50, we obtain values of — 0.15 compatible with the WMAP result 
(tt = 0.17 ± 0.06 for the combined fit). The empirical 2a limit, tt > 0.05, implies Az > 1. 
In the figures, we will focus on the Az = 1 case since it is close to the minimum ionizing 




(11) 
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photon requirement B — 1 (obtained in the absence of recombinations), while remaining 
marginally compatible with the WMAP results. 

The numbers given in Table 1 can be used directly to obtain the luminosity of the 
sources of ionizing photons on the assumption that they are a population of identical objects 
with a specified surface number density. Indeed, the numbers in the table represent the 
apparent magnitudes of individual objects if their surface number density is 1 arcmin^^. A 
simple rescaling yields the apparent magnitude m for different assumed surface densities J\f 
namely 

log^^ ^0.4:{m- hab), (12) 

where /iab is the mean surface brightness given in Table 1. This relationship for Z2 — 6, 
Az — 1, and an observer wavelength A ~ 1400(1 + Z2) A is shown in the left panel of 
Figure 1 as the solid line labeled (1, 1) (this notation will be explained in the next Section). 
Models characterized by different values of mean surface brightness are obtained by simple 
translations with respect to this reference model. 

Note that in the idealized case of identical sources, an entire population of such sources 
is represented by a single point in diagrams such as Figure 1. Those lying on the indicated 
line in the diagram would be sufficient to ionize the IGM. Observations that fail to detect 
sources at some surface density-luminosity combination along the line would not rule out 
the model but instead could indicate that the sources lie at some other point along the line. 
The shaded area below the fine indicates the area where sources would not produce enough 
ionizing photons for reionization. 

The results for a different effective temperature arc illustrated in the right panel of 
Figure 1, which is analogous to the left panel but refers to sources with = 5 x 10^ K. 

2.5. Mean metallicity at reionization 

It is interesting to compute the mean metallicity to which the universe is enriched by 
stellar populations with the minimum mean surface brightness required to ionize the IGM. 
Here we define metallicity as Z = Qz/^b, where Qz is the fractional comoving density in 
metals (in both galaxies and the IGM). Massive stars with effective temperatures in the 
range 3 x 10^ to 10^ K produce 10^^ Lyman continuum photons per unit of stellar mass. 
Assuming a typical yield in metals of 0.2 (Woosley and Weaver 1995), we find that the 
massive stars responsible for reionization enrich the universe to a minimum metallicity of 
Zjjiin ~ 1.2 X 10~^ Zq. This result does not change by more than a factor 3 for a wide range 
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of stellar properties. As an example, by considering only supermassive Pop III stars, with 
masses between 100 and 1000 Mq with the yields as given by Hcgcr and Woosley (2002), 
and a power- law mass function with a slope between —1 and —2, we find that the mean 
mass-weighted yield is ~ 0.2 ± 0.05 (neglecting the metals trapped in black hole remnants). 
Stars in this mass range will produce 7.6 x 10^^ ph Mq^ (Tumlinson ct al. 2001), so 
that the minimum surface brightness model enriches the universe to a mean metallicity 
'Z'min ~ 1-7 X 10~^Zq. For a given IMF, any reduction of the output of ionizing photons 
can be regarded as an inefficiency in the energetics of reionization and will increase the 
fraction of hydrogen that needs to undergo nuclear processing and therefore will increase the 
mean metallicity. Thus, the maximum efficiency, minimum surface brightness model also 
corresponds to the minimum metallicity case. 

Estimates of the mean density of heavy elements in the universe at the end of reionization 
{z ~ 6) place additional constraints on the mean surface brightness of the stellar component 
of reionization sources. It is difficult to determine flz precisely since some metals will be in 
stars, some in the interstellar media of galaxies, and some in the IGM. But we can bracket 
the likely range of metal densities as follows. Songaila (2001) has estimated a mean metal 
density in the IGM of ^z,igm ~ 3 x 10^^, corresponding to Z ~ 3 x IO^^Zq, at 2; ~ 5 from 
the statistics of CIV and SilV absorption-line systems. This provides a firm lower limit to 
the mean metal density shortly after the end of reionization. This limit is about twice the 
metallicity of the minimum surface brightness model Zmin (discussed above). 

The present mean metallicity in galaxies is close to solar (Kulkarni and Fall 2002) and 
the present mean density of the stellar and interstellar components of galaxies is flgai ~ 
4 X 10~^ (Fukugita et al. 1998). Together these imply flz,gat ^ 10~^ at 2; = 0. A variety 
of observations and models indicate that the fraction of mass in galaxies ftgai was smaller 
at 2; ~ 6 than at 2; = by a factor of 10 or more, and that the mean metallicity in 
galaxies was also smaller by a similar factor (Pei et al. 1999). From these considerations 
we infer a generous but non-rigorous upper limit on the mean metal density of flz < 10~^, 
corresponding to Zmax = ^z/^b < O.OIZ0 at 2; ~ 6. This value is ~ 60 times higher than 
the metallicity of the minimum surface brightness model. 

If metal-producing sources are surrounded by local HII regions, their luminosity at 
A = 1400 A can be enhanced by up to a factor 3 due to the nebular continuum contribution 
(see Section 3.1). Including this factor, we use the limit to metal production as an upper 
limit to the mean surface brightness of reionization sources. This is obtained for Z2 = Q and 
Az = 1. Reionization at higher redshifts or occurring over a more extended time interval, 
produces the same (or higher) metallicity for fainter sources. Thus, the regions under the 
global metaUicity constraint lines in Figure 1 define necessary but not sufficient conditions 
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to guarantee that metals are not overproduced. All reionization models compatible with 
the global metallicity constraint must lie in the non-shaded region in Figure 1. Our results 
are broadly compatible with those of Miralda-Escude and Rees (1998) once the different 
assumptions are taken into account. 

3. More realistic models 

3.1. Escape of ionizing photons 

It is likely that some Lyman continuum photons are absorbed within the reionization 
sources themselves. Even in an environment with primordial composition, and thus free of 
dust, a lower escape fraction fc can result if the individual sources are surrounded by an 
envelope of neutral hydrogen left by incomplete star formation in protogalaxies. In such 
clouds, recombinations would occur at much higher rates; these objects would produce local 
HII regions for at least part of their lifetimes and thus could display hydrogen and helium 
emission lines in their spectra (Bromm et al. 2002; Haiman 2002). We do not know the 
best value of the mean escape fraction. It only takes a column density of 1.7 x 10^'' hydrogen 
atoms cm~^ to produce unit optical depth at the Lyman limit (A = 912 A). Moreover, 
some amounts of dust may well be present because of enrichment by previous generations 
of stars. Steidel et al. (2001) derive from their observations, corrected for intervening IGM 
absorption, a reduction of the Lyman continuum with respect to the continuum at 1500 A 
by ~ 50 %. They estimate that dust absorption in the far-UV be a factor of ~ 5. Thus, the 
total combined effect of gas and dust gives fc ~ 0.1. The Steidel et al. (2001) observations 
also suggest that metal-free, dust-free galaxies might have fc ~ 0.5. However, the properties 
of the reionization sources at z > 6, including the Lyman continuum escape fraction, may be 
different from those of the Lyman-brcak galaxies at 2; ~ 3. Leitherer et al. (1995) argue for 
escape fractions below 15% in bright starburst galaxies in the local universe. However, these 
observational studies measure the escape fraction near 912 A while the hydrogen ionization 
cross section decreases rather rapidly with decreasing wavelength so that even with optical 
depth T = 2 — 3 at 912 A one could have a Lyman continuum escape fraction greater than 
20 %, especially for supermassive stars that are hotter than local O stars (Mezger 1978; 
Panagia et al. 2003). There have also been a few attempts to estimate the escape fraction 
theoretically based on idealized assumptions on the structure of the sources. The derived 
values of the escape fraction span a relatively wide range roughly centered around fc — 0.1 
(Dove and Shull 1994; Dove et al. 2000; Ricotti and Shull 2000). In the absence of a 
definite value, we will consider two possible escape fractions: fc = 0.5 and fc = 0.1. 

An escape fraction fc<i increases the required mean surface brightness relative to the 
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minimum value. This in turn increases the required value of S by a factor f~^. Moreover, 
the fraction of Lyman-continuum radiation that does not escape generates a nebular contin- 
uum (mostly two photon emission) in the rest-frame non-ionizing UV, which dominates the 
observed continuum flux at A = 1400 A (Panagia et al. 2003). Combining the two effects, 
we find that the actual continuum flux F at A = 1400 A is higher than the value for the 
complete escape of ionizing photons by the factor 

EifA ^ (13) 

F(l) - fc ' ^'^^ 

where the coefficients 3 and 2 account for the fact that the nebular continuum at 1400 A for 
population III HII regions is roughly twice as high as the stellar continuum. 



3.2. Clumpy IGM 

In addition to inhomogeneities near and within the sources, there will also be density 
fluctuations on larger scales. Such inhomogeneities will increase the recombination rate and 
act as sinks of ionizing UV photons. This effect can be included by changing the value of the 
dumpiness factor C in equation (9). We have considered the value of C = 30 (Gnedin and 
Ostriker 1997; Madau et al. 1999) and recomputed the mean surface brightness and the 
Thomson optical depth tt for the same models that we considered for C = 1. The resulting 
values are also given in Table 1. Based on the simulations by Gnedin and Ostriker (1997), 
the value of C = 30 appears to be a reasonable choice for 2; ~ 6 — 7. At higher rcdshifts, 
smaller values would be expected. Thus, it is likely that the effective value of C is bracketed 
by the two values (C = 1 and 30) we consider here. 

As a result of gravitational instability, more reionization sources will form and the 
IGM will become more clumpy as time increases. Thus, the increase in the ionization and 
recombination rates will at least partially offset each other. This argument suggests that our 
results, based on the idealization that fic and C are constant, will be more accurate than 
they might at first appear. 

An interpolation formula connecting the parameters B and C is given in the Appendix. 
We find that the value of C does not affect the optical depth tt significantly. In contrast, 
the surface brightness increases by a factor of 2-10 (i.e. 0.8-2.5 magnitudes) for the models 
with C = 30, as shown in Table 1. 



3.3. Two-parameter models 



In the following, we denote models with the notation (/c, C) to indicate the adopted 
values of the Lyman-continuum escape fraction and of the dumpiness parameter. The (1,1) 
label indicates the minimum surface brightness model as described in the previous section. 
The B factor of equation (5) is now a function of fc and C. Any model with C — 1 located 
below the metallicity constraint in Figure 1 is allowed. Note that the maximum-metallicity 
constraint was derived in the limit of very small fc so that the nebular continuum can enhance 
the source luminosity by up to a factor of three. Models with f,- — 1 and C > 1 would 
have little or no contribution from nebular continuum but would be intrinsically bright and 
therefore produce more metals. Thus, it is possible that some models with C > 1 and /c — 1 
lie below the shaded area but still violate the maximum metallicity constraint. In Figure 1, 
we plot as thin dotted lines the luminosity-surface density relations for identical reionization 
sources from models with (/c, C) = (0.5,1), and (0.1,30). The lines were computed for 
Z2 — 6 and Az — 1 but are nearly identical for all 2^2 < 10. The mean surface brightness of 
the reionization sources changes with Az by less than a factor 4 for all C < 30 (see Table 
1). For the sake of completeness, we have explored a wider range of parameters than it 
is probably realistic. Note that both of the (0.5,1) and (0.1,30) models are allowed by the 
maximum metallicity constraint. 



3.4. Luminosity function of ionizing sources 

So far, we have considered reionization sources with identical luminosity, but it is more 
likely that they are characterized by a broad luminosity function. On the basis of both theo- 
retical considerations and experience at lower redshifts, it is reasonable to adopt a Schechter 
(1976) differential luminosity function, $(-L) oc exp (— L/L*). Here, we adopt the pa- 
rameters derived for Lyman-break galaxies at 2; ~ 3 — 4, i.e. a = 1.6 and M^^i4oo = —21.2 
as our reference luminosity function (Steidel et al. 1999; Yan et al. 2002) and explore the 
effects of variations in both M^.^i4oo and a. The top left panel of Figure 2 was computed for 
the reference cumulative luminosity function in the case of 2^2 = 6 and Az — 1. The bottom 
solid line refers to the minimum surface brightness (1,1) model. The area beneath this line 
would not reionize the IGM. The top sohd line shows the metaUicity constraints and the 
shaded area above this line shows the region where - with the adopted luminosity function 
- metals would be overproduced. The (0.5,1) and (0.1,30) models are shown as dotted lines. 
The remaining panels of Figure 2 explore the effects of varying 1400- The top right panel 
is for M^, 1400 = -23.2, i.e. two magnitudes of luminosity evolution compared to z = 3 Ly- 
man break galaxies. The bottom left panel is for 1400 = -17.5 corresponding to a star 
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cluster with IO^Mq in massive stars forming over a period of ~ 1 Myr. Finally, the bottom 
right panel is for 1400 = —13.5 which represents 7.7 magnitudes of luminosity evolution 

compared to the z = 3 Lyman break galaxies and corresponds to the case of small dwarf 
galaxies or proto-globular clusters. The non-shaded area in the figures indicates the allowed 
region for each luminosity function. The curves would shift by less than 0.1 magnitudes for 
any value of Z2 < 10. 

In contrast to those in Figure 1, the lines in Figure 2 represent luminosity functions that 
need to be populated everywhere in order to produce the required ionizing flux. Thus, a sur- 
vey ruling out sources at one point along a curve would rule out that particular combination 
of (/c, C) model and luminosity function. 

In Figure 3, we explore the effects of changing the slope of the luminosity function. The 
top left panel is identical to that of Figure 2 and is repeated here for convenience. The top 
right panel is for a = 1.1, similar to the slope of the z = luminosity function. The bottom 
left panel is for a — 1.9. The bottom right panel of Figure 3 shows the forbidden region 
(double shading) and the preferred region (no shading) for a range of models. In this last 
panel, we have considered models with 1400 ranging from -21.2 to -16.2 and luminosity 
function slopes ranging from a — 1.1 to 1.9. The lower shaded area represents the region 
of underproduction of ionizing photons common to all these luminosity functions. Similarly, 
the upper shaded area corresponds to the common region of overproduction of metals. The 
non-shaded area is allowed for all the models within the adopted range of parameters. Thus, 
this area should be favored in searches for the reionization sources. 

The analogs of Figures 2 and 3 for models other than those with Z2 = Q and Az = 1 need 
in general to be recomputed by properly renormalizing the luminosity function. However, 
these curves can be obtained approximately by translating vertically the curves for the (1,1) 
model with z = 6 and Az = 1 by an amount AlogM = —OAAfiAB, where Afi^B is the 
difference between the values of mean surface brightness of the models considered. This 
procedure is only approximate because it does not properly average the value of M*_i4oo over 
the redshift interval Az being considered. 

3.5. Detecting Lya from ionizing sources 

The fraction 1 — /c of the ionizing continuum that does not escape generates local HII 
regions around the reionization sources, making them potentially visible in the H recombi- 
nation lines and, more faintly, in the Hel lines. Prior to full reionization, when the universe 
consisted of sources embedded in their own HII regions, Lya photons were able to escape if 
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these HIT regions were big enough that their outer boundaries were redshifted out of reso- 
nance (Madau and Rees 2000; Haiman 2002). In the absence of dust, the fraction of Lya 
photons that escape from the HII region is determined solely by its size and hence by the 
production rate of ionizing photons of its central source Nc (Loeb and Rybicki 1999). In 
the more realistic case, where some dust is present, an even smaller fraction of Lya photons 
would escape (Panagia and Ranieri 1973; Bonilha et al. 1979; Chariot and Fall 1993). 
Neglecting absorption by the neutral medium, the Lyct line intensity can be computed as fol- 
lows. Each recombination within the HII region has a 69% probability of producing one decay 
leading to Lya emission. Under conditions of very low metallicity {Z < 10'^ Zq), the Lya 
production is increased by ~50 % due to coUisional excitation (Panagia et al. 2003). The 
number density of ionizing photons absorbed within the local HII region is obtained from the 
required number density of photons needed to reionize the universe, as (/~^ — l)nc, where 
fc is again the escape fraction of Lyman continuum photons. Following Miralda-Escude and 
Rees (1998) - and in analogy to equation (5) - we find that the rate of production of Lya 
photons per unit solid angle is: 



The Lya surface brightness is the product of dNa/dfl and the mean observer frame photon 
energy hPa- For Az — 1 and Z2 — 6, this surface brightness is {f~^ — 1) x 1-276 x 10~^^ 
erg cm~^ s~^ arcmin"^. This value is a strong function of C, decreases with increasing Z2, 
and depends only slightly on the duration of reionization. We now define fa as the escape 
fraction of Lya photons. The average surface brightness in Lya, J^, of the reionization 
sources corresponding to the {fc, C) model is then expressed as 



J, = ^^(/c"' - i)fahPanH. (15) 

We have computed fa following Miralda-Escude (1998) and Madau and Rees (2000), as- 
suming a line width of 150 km s^^. In the left panel of Figure 4, we show the locus of 
identical Lya sources corresponding to {fc,C) = (0.5,1) and (0.1,30) and for Z2 = 6 and 
Az — 1. We also show the metallicity constraint as the upper shaded area. In the case 
of Lya sources, there is no lower shaded area, indicating insufficient production of ionizing 
photons, as sources with fc close to unity have arbitrarily low Lya flux but are still capable 
of ionizing the IGM. The right panel of Figure 4 shows two luminosity functions for the Lya 
sources. The solid line represents the same luminosity function as the z = 3 Lyman-break 
galaxies converted to Lya luminosity for the (0.5,1) model. The dashed line represents a 
luminosity function with a — —1.1 and 1400 = —17.5. 
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In addition to the case of an intrinsic width of 150 km s~^, we have also considered cases 
with 75 to 300 km s"-'^. As expected, the wider the hne, the more Lja radiation escapes. 
For example, for sources at the faint end of Figure 4, the fraction of escaping Lya photons 
is ~ 30 % higher for v — 200 km s~^ and ~ 30% lower for v = 100 km s~^. 

4. Present and Future Surveys 

Our models can be used to plan and interpret searches for reionization sources with 

present and future telescopes. A possible approach to identifying the reionization sources 
is the Lyman-break technique, which consists of broad-band photometry with three filters, 
two to measure the continuum longwards of Lya and a third filter to measure (or place 
upper limits on) a much lower fiux just below Lja. Note that for z <4, the technique 
places the shortest wavelength filters across the Lyman-continuum break, while at z ^5 the 
combined effect of the Lya forest shifts the wavelength of the break to Lya. In Figures 1, 
2 and 3, we identify the region in the surface density vs luminosity plane probed by the 
GOODS ACS survey (Dickinson and Giavahsco 2003), the HDF (Thompson et al. 1999) 
and HDFS (WiUiams et al. 2000) NICMOS fields, the Hubble Ultra Deep Field (UDF) 
and a hypothetical survey with the James Webb Space Telescope (JWST) with 100 hr of 
integration per filter. The limiting magnitudes are those corresponding to a lOa detection 
in the prescribed integration time. The limiting surface densities of these surveys has been 
determined by requiring that at least three sources be detected at or above lOcx significance 
level in the survey area. Thus, the limiting surface density is given by N'nm = 3/ (survey area). 
These two limits define the quadrant in the logAf — ttiab plane where a particular survey 
is sensitive to reionization sources, i.e. higher Af and brighter magnitudes. Note that the 
curves in Figure 1 again apply to the case of identical sources. Thus, a survey not finding 
sources at its limiting depth and surface density does not rule out a particular model but 
only a particular luminosity-surface density combination. In contrast, in Figures 2 and 3 
where broad luminosity functions arc plotted, all sources along the particular curve being 
considered arc needed to rcionizc the IGM. Thus, a survey not finding sources at its limiting 
luminosity and surface density can rule out that model-luminosity function combination. 

It is clear that GOODS/ACS observations cannot reach the faintest predictions shown 
in the figures. Depending on the actual luminosity function parameters a and M* i4oo, 
GOODS/ACS may be more or less successful than the UDF in identifying z = 6 reionization 
sources. In general, ACS observations could identify reionization sources with the F850LP 
filter if reionization occurred at a low redshift (z~6) and if the reionization sources had a 
luminosity function similar to that of Lyman-break galaxies. However, even for Z2 — 6, one 
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can detect all reionization sources as F775W dropouts only if ^ 0.5. For larger values 
of Az or higher values of the reionization redshift, many sources would drop out of F850LP 
and could be identified only by a deep follow-up study, e.g., with the IR channel of WFC3 
(Stiavelli et al. 2000). While the HDF/HDFS NICMOS observations could be used to 
probe higher redshifts, they lack the combination of area and sensitivity necessary to detect 
reionization sources as faint as the minimum surface brightness ones. In contrast, JWST 
has both the wavelength coverage and the sensitivity to detect minimum surface brightness 
reionization sources for all the models discussed here, provided that M*_i4oo < —14.2; see, 
for example, the bottom panels in Figure 2. 

Another approach to search for the reionization sources before the advent of JWST 
is to use gravitational lensing by foreground clusters of galaxies to increase the sensitivity 
of the surveys (Ellis et al. 2001; Hu et al. 2002). In Figure 5, the hne identified as 
"Lens" represents a hypothetical survey with 100 orbits per filter with the ACS. The field is 
centered on one cluster of galaxies at z = 0.5, idealized as a singular isothermal sphere with 
a one-dimensional velocity dispersion of 1000 km s^^, corresponding to Tx — 10 keV. The 
faint-end of the line represents an amplification by a factor 10. The slope of this line in the 
intrinsic density of sources (in the source plane) vs. flux plane is -2. This slope corresponds 
to a constant observed surface brightness since the apparent (i.e. in the lens-plane) density 
of sources is lower by a factor identical to the gravitational amphfication. The luminosity 
function assumed for this figure has a slope a — —1.6 and a knee at M*^i4oo = —17.5. Clearly, 
these sources would be too faint to detect in a survey like the UDF but would be detectable 
in a 400 orbits survey with ACS pointed toward a suitably chosen cluster of galaxies acting 
as a lens. Such a survey could probe an interesting area of parameter space, especially if 
foUowed-up by WFC3 near-IR imaging to extend its sensitivity to higher redshifts. 

An alternative approach to the Lyman-break technique is the direct detection of Lya 
emission. As pointed out by Haiman (2002) and others and discussed in Section 3, Lya 
emission can be detected, with less sensitivity, even from sources at redshifts before reion- 
ization is completed (i.e., z > Z2) as long as either the line is intrinsically wide or the source 
is bright enough to generate its own "proximity effect" by ionizing neutral hydrogen in its 
neighborhood. Figure 4 shows the detection limits for a survey based on the narrow-band 
excess technique (Rhoads and Malhotra 2001; Hu et al. 2002), which consists in obtaining 
narrow-band exposures to detect line emission and broad-band images or adjacent narrow- 
band exposures to measure the continuum level. We assume that the survey would be 
carried out at 2; > 6 using a large field-of-view camera on a 4-meter ground-based telescope. 
In the case illustrated in Figure 4, the search is assumed to cover ~ 0.6 square degrees with 
AA/A ^ 0.01 (labelled LALA in the Figure). Since Figure 4 shows the density of all Lya 
sources, the effective area explored by a narrow-band excess survey is reduced by a factor 
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proportional to the spanned by the filter. The obhque marker identified as Lens in Figure 
4 represents the detection hmits of a 100-orbit ACS grism survey exploiting gravitational 
lensing. In addition, we plot in Figure 4 the upper limits (arrows) and the single detection 
(solid vertical bar) reported by Hu ct al. (2002). Taken at face value, these upper limits 
would rule out, for the present model, values of M* brighter than ~ —17.5, thus suggesting 
that only dwarf galaxies exist before reionization. The estimated density of objects like the 
one detected by Hu et al. (2002) at 2; ~ 6.56, by gravitational lensing amplification, would 
be sufficient to reionize the universe if these objects had hot ionizing continua with a high 
escape fraction, 0.1 < /c < 0.5. Clearly, before strong conclusion can be drawn, we need to 
identify and study more sources at 2; > 6. 



5. Conclusions 

We have developed a method to predict the observability of cosmological reionization 
sources based on the fact that reionization requires a known comoving density of ionizing 
photons. This number needs to be at least one per atom, thus providing a model-insensitive 
minimum mean surface brightness (in fiux, maximum in magnitudes) for the reionization 
sources. This minimum mean surface brightness of reionization sources is general and robust. 
Our particular choice of cosmological parameters does not affect our results by more than a 
few tenths of a magnitude. We find that if reionization is caused by UV-efficient, minimum 
surface brightness sources, the non-ionizing continuum emission from reionization sources 
will be difficult to detect before the advent of JWST. However, the actual mean surface 
brightness may be many times higher than this ideal limit if some UV photons are absorbed 
within the sources themselves and/or if some ionizations are offset by recombinations in the 
clumpy intergalactic medium. In this more realistic case, HST has a reasonable chance of 
detecting the reionization sources in very deep images, especially if the sources are observed 
behind clusters of galaxies, so that they are magnified gravitationally. If instead the sources 
of reionization were not extremely hot population III stars but cooler population II stars or 
AGNs, they would be brighter by 1-2 magnitudes and thus they would be easier to detect. 

The predictions for Lya emission intensities are less secure than our continuum fluxes 
because they depend on geometry and line width and could be diminished by absorption by 
dust, which we have ignored here for simplicity. However, if the Lya intensities are close to 
our estimates, searches based on narrow-band excess techniques or slitless grisms would be 
promising and might lead to the detection of the reionization sources within this decade. 

We wish to thank Tommaso Treu and Mike Gladders for helpful discussions. Comments 
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from an anonymous referee contributed to significantly improving the presentation of our 
results. This work is partially supported by NASA Grant NAGS- 12458. 

A. An interpolation formula for B{C) 

As discussed in Sections 2 and 3 (see also Table 1) the factor B, which provides a measure 
of the extra ionizations needed to compensate for recombinations and for competing He 
ionizations, is a sensitive function of both the dumpiness factor C and the redshift interval 
Az over which the universe rcionization takes place. This is because both a high-rcdshift 
and a large C imply higher densities and, therefore, an overall higher recombination rate 
that requires a larger number of ionizing photons. 

On the other hand, the direct dependence of B on the reionization redshift for any fixed 
Az is much weaker, because when the reionization process is assumed to take place over a 
fixed Az interval, for each redshift there is partial compensation between a higher density 
oc {1 + z)^ and a shorter duration of the phenomenon. At oc (1 + z)~^/'^. 

These facts are clearly illustrated in Figure 6, in which the behaviour of B as a function 
of C is shown for Z2 — 6 and several values of (left-hand panel), and for Az — 1 and 
several values of the reionization redshift (right-hand panel) . 

Given the importance of the parameter B in cosmological applications, it is convenient 

to devise a simple formula that allows a fairly accurate and quick calculation of B for any 
value of z, Az, and C, while retaining the correct functional dependences in the asymptotic 
limits. After a little experimentation, we find the following approximate formulae: 

B = 1.09[1 + {^CfT^^, (Al) 

^ = 0A/[i + {6/Az){6/z2n m 

p = 0.30 + 0.6(V6)(Az/6). (A3) 

These formulae correctly tend to a constant value of B=1.09 in the limit of small values 
of C and to B directly proportional to C for high C values. Equation (Al) gives an excellent 
fit (better than 6%) for Z2 = 6 and all values of Az as shown in Figure 6 (left-hand panel), 
as well as for cases of Az — 1 and any value of Z2 (right-hand panel). For all other cases, the 



-19- 



accuracy is still within about 12%, which is adequate for most applications. The fit can be 
optimized for specific needs by adjusting the coefficient of the z-it^z term in equation (A3). 
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Fig. 1. — Surface density vs apparent AB magnitude for identical reionization sources at 
Z2 = 6, Az = 1. A line represents the locus of identical sources with a given mean sur- 
face brightness. The bottom solid line represents the minimum surface brightness model, 
(/c, C) = (1,1), while the upper solid line represents the global metallicity constraint 
Z < 0.01^0 at z = Z2. The thin dotted lines represent the (0.5, 1) and (0.1, 30) mod- 
els which are characterized by a lower escape fraction of ionizing photons (see text). The 
non-shaded area is the only one accessible to reionization sources that do not overproduce 
metals. The L-shaped markers dehmit the areas probed by the GOODS/ACS survey, the 
HDF/HDFS NICMOS fields, the Ultra Deep Field (UDF) and by a hypothetical ultra-dccp 
survey with JWST. The left panel refers to pop 111 sources with effective temperature of 10^ 
K. The right panel to reionization sources with effective temperature 5 x 10^ K. 
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Fig. 2. — Surface density vs apparent AB magnitude of reionization sources with luminosity 
functions with different knees. The lower solid line represents the minimum surface brightness 
model, (1,1), for = 6 and Az = 1. The top thin solid line represents the global metallicity 
constraint. The thin dotted lines give the luminosity function for the (0.5,1) and the (0.1, 
30) models (see text). The non-shaded area is the only one accessible to reionization sources 
with the specified luminosity function. The top left panel refers to the z — 3 Lyman-break 
luminosity function, the remaining panels explore the effect of changing the knee M*^i4oo. 
The L-shaped markers delimit the areas probed by the GOODS/ACS survey, the HDF and 
HDFS NICMOS fields, the UDF, and an ultra-deep survey with JWST. If the value of M.^mqo 
were as faint as -13.5, the reionization sources would be hard to detect even with JWST. 
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Fig. 3. — Surface density vs apparent AB magnitude of reionization sources with luminosity 
functions with different faint-end slopes. The lower solid line represents the minimum surface 
brightness model, (1,1), with Z2 = 6 and Az = 1. The top thin solid line represents the 
global metallicity constraint. The thin dotted lines give the luminosity function for the 
(0.5,1) and the (0.1, 30) models (see text). The non-shaded area is the only one accessible 
to reionization sources with the specified luminosity function. The top left panel refers to 
the z — 3 Lyman-break luminosity function, the top right and bottom left panels explore 
the effect of changing the faint-end slope a. The bottom right panel shows the forbidden 
region (doubly shaded) and the preferred region (non-shaded), corresponding to luminosity 
functions with a ranging between 1.1 and 1.6 and M^, 1400 ranging between -16.2 and -21.2. 
The L-shaped markers delimit the areas probed by the GOODS/ACS survey, the HDF and 
HDFS NICMOS fields, the UDF, and an ultra-deep survey with JWST. 
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Fig. 4. — Surface density vs Lya line flux of reionization sources. The thin dotted hues in 
the left panel indicate for two different (/c, C) models the locus of identical sources for Z2 = Q 
and Az = 1 (see text). The top solid line identifies the global metallicity constraint. The 
luminosity function adopted in the figure has a = 1.6 and M^, = —16.2. The right panels 
illustrates two different luminosity functions. The solid line refer to a luminosity function 
identical to that of 2; = 3 Lyman break galaxies, while the dashed line refers to a luminosity 
function with M*^i4oo = —17.5 and a slope identical to the local slope, o; = 1.1. In both 
panels, the L-shaped marker identifies the area probed by a narrow-band excess survey at 
z > 6. The oblique marker labeled Lens represents a hypothetical 100-orbit survey with 
the ACS grism on a cluster of galaxies to exploit gravitational amplification. The solid bar 
represents the density estimated from the detection at z = 6.56 by Hu et al. (2002) while 
the two points with down-pointing arrows represent their upper limits. 
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Fig. 5. — Surface density vs apparent AB magnitude of reionization sources with a broad 
luminosity function. The lower solid line represents the minimum surface brightness, (1,1), 
model with Z2 = 6 and Az = 1. The top thin solid line represents the global metallicity 
constraint. The thin dotted lines give the predictions for the (0.5,1) and the (0.1, 30) models 
(not labeled, see text). The adopted luminosity function has faint-end slope a — —1.6 
and knee 1400 = —17.5. The non-shaded area is the only one available to reionization 
sources with this luminosity function. The L-shaped markers delimit the areas probed by the 
GOODS/ACS survey, the HDF and HDFS NICMOS fields, and an ultra-deep survey with 
JWST. The oblique marker labeled Lens shows the area probed by a 400-orbits survey with 
ACS pointed toward a cluster of galaxies as gravitational lens. If the value of 1400 is as 
faint as -17.5, ACS will be able to detect the reionization sources only if they are relatively 
inefficient and only if the surveys exploits gravitational amplification. 
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Fig. 6. — The left panel shows the mean number of photons per ionization S as a function 
of the IGM dumpiness factor C for — & and a variety of values of Az. The solid squares 
were obtained by direct integration of equation (8), while the lines have been obtained with 
equation (Al). The right panel shows S as a function of C for a given value of A2; and 
different values of Z2. The variations are more modest because, for a given A2;, the shorter 
time interval at increasing Z2 is compensated by the higher recombination rate due to the 
higher density. 
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Table 1. Mean surface brightness for reionization for H+He IGM and = 10^ K 





Az 


C 


B 


TT 


Ml400{l+z2) 
(AB mag arcmin^'^) 


A'l40()(l+z) 
(AB mag arcmin"^) 


^Lya(/c =0.5) 

(erg cm s^^ xlO^®) 


6.2 


1 


1 


1.115 


0.046 


28.773 


28.874 


1.276 


6.2 


1 


30 


2.267 


0.047 


28.003 


28.103 


2.594 


6.2 


3 


1 


1.158 


0.054 


29.518 


28.929 


1.208 


6.2 


3 


30 


4.998 


0.062 


27.930 


27.341 


5.216 


6.2 


10 


1 


1.260 
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30.238 


29.240 


1.119 


6.2 


10 


30 


9.042 


0.097 


28.099 


27.101 


8.031 


6.2 


20 


1 


1.335 


0.105 


30.439 


29.507 


1.088 


6.2 


20 


30 


10.732 


0.128 


28.176 


27.245 


8.748 


6.2 


30 


1 


1.373 
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30.063 


29.716 


1.078 


6.2 


30 


30 


11.373 


0.150 


27.768 


27.420 


8.930 


6.2 


40 


1 


1.397 


0.140 


30.349 


29.941 


1.075 


6.2 


40 


30 


11.694 


0.166 


28.043 


27.635 


8.999 


6.2 


50 


1 


1.413 


0.153 


30.131 
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